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Abstract. The motion of free nearby test particles relative to a stable equatorial 
circular geodesic orbit about a Kerr source is investigated. It is shown that the 
nonlinear generalized Jacobi equation can be transformed in this case to an autonomous 
form. Tidal dynamics beyond the critical speed c/\/2 is studied. We show, in 
particular, that a free test particle vertically launched from the circular orbit parallel 
or antiparallel to the Kerr rotation axis is tidally accelerated if its initial relative speed 
exceeds c/a/2- Possible applications of our results to high-energy astrophysics are 
briefly mentioned. 
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1. Introduction 

The purpose of this paper is to study some aspects of the relative motion of free test 
particles in the gravitational field of a rotating astronomical source. Imagine, for the 
sake of concreteness, a free mass m in orbit about the central source of mass M ^ m 
and angular momentum J, where m, M and J are constants. Let r be the proper time 
along the orbit of mass m and A^j-^-j be a local orthonormal tetrad system that is parallel 
transported along this path. Here A'^j.q^ = dx^/dr is the four- velocity vector of m and 
^^(j^, = I5 2, 3, are unit spacehke gyro directions that determine the local spatial frame 
of m. We employ units such that c = G = 1 throughout; moreover, the signature of 
the metric is +2 in our convention. Suppose that at r = a probe is launched from 
m and for some time r > the motion of the probe is closely monitored by observers 
comoving with m. To study relative motion in an invariant manner within the context of 
general relativity, the quasi-inertial Fermi normal coordinate system pi is indispensable 
as it most closely corresponds to actual observations. Let X" = (T, X) be the Fermi 
coordinates of the probe at an event P, then there exists a unique spacehke geodesic 
of length £ that connects P orthogonally to the path of m at an event O with proper 
time T such that T = t and Xi = in^^X'^^.y where = {dx^/di)o is the unit vector at 
O that is tangent to the spacelike geodesic at O and orthogonal to the worldline of m. 
The metric in Fermi coordinates is given by 

^700 = - 1 - ^Ro^OJX'X^ + ■ ■ ■ , (1) 
2 

9oi = — 2 RojikX^X + ■ ■ ■ , (2) 

9ij = ^ij ~ 2 ^ RikjlX^ + • • • (3) 

to second order in the distance away from m that permanently lies at the spatial origin 
of the Fermi coordinates by construction and has Fermi coordinates (r, 0). Here ^Rap-yS 
is defined by 

which is the projection of the Riemann curvature tensor on the tetrad frame A^j.^^ along 
the reference worldline. The theoretical equation of motion of the probe relative to mass 
m can be described by the reduced geodesic equation [2] 

^ + (r-,-r»,n— — = 0. (6) 

where V = dl^/dT is the Fermi velocity of the probe relative to m. The four- velocity 
of the probe is given by = r(l, V), where F can be determined from 

-^=goo + 2goX' + gijV'V^ < 0, (6) 

since f/^ is a timelike unit vector. The Fermi coordinates are geodesic coordinates based 
on an orthonormal tetrad frame A'^^^^(r) that is parallel propagated along the orbit of 
m. They are admissible within a cylindrical region of radius TZa along the worldline 
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of m; that is, |X| < TZa, where TZaiT) is a certain radius of curvature of spacetime. 
The nature of TZa{T) has been discussed in detail in our recent paper on exphcit Fermi 
coordinates 

For the sake of simphcity, we hmit our considerations to |X| ^ TZa- A general 
feature of the reduced geodesic © is that it contains a critical speed given by l/\/2, 
especially in the case of one-dimensional motion. The notion of critical speed in 
gravitational motion has been recently reviewed in ref . |S| . We are particularly interested 
in the role that the critical speed Vc = l/\/2 plays in the motion of the probe; therefore 
the probe can be launched from m at any initial speed Vq < 1, but — for simplicity — we 
concentrate on its motion only when it is relatively close to the reference particle m. 
Equation (0) to first order in |X| is the generalized Jacobi equation |2] 



dT^ 

+l{3^RokjoV'X' + ^R^kjiVv' + ^RokjiV'V'V')X^ = 0, (7) 
which describes the timelike motion of the probe when 

-L = 1 _ 1/2 + ^R,,,^X'X^ + ^^i?o,,fcX^V^*X^- + ^^R,k,iV'X''V^X' > 0. (8) 

For |V| ^ 1, we can drop all velocity-dependent terms in equations ((Tj) and (jH)); then, 
the inequality in (jH} is satisfied and ((Zj) reduces to the Jacobi equation 

^ + ^Ro^ojX^ = 0. (9) 

Suppose that the worldline of the probe can be determined in the Fermi coordinate 
system on the basis of equations (0) and ©; then, it is in principle possible to find the 
representation of this worldline in the background coordinate system using the explicit 
coordinate transformation between the two systems of coordinates. In practice, however, 
it is not possible in general to find explicit expressions for Fermi coordinates in terms of 
the background coordinates and therefore one has to resort to approximation schemes; 
see ref. j2] for a general discussion of this problem. Nevertheless, it is important to 
recognize that once X(T) is determined from the generalized Jacobi equation (or the 
Jacobi equation), the worldline of the probe is completely characterized. That is, if 
in the background coordinate system the worldline of m is given by x^{t), then the 
worldline of the probe is xfi^r), where 



x;^(r)=x'^(r)+X^(r)A'^(,)(r) (10) 

and we have replaced the Fermi temporal coordinate T by r in the argument of 
X. This relation follows from the definition of Fermi coordinates and the fact that 
x!^{t) — x^{t) ^ in^ for i <^ TZa, since the (generalized) Jacobi equation is linear in 
the distance away from the reference trajectory. Alternatively, x^{t) may be viewed as 
a perturbed geodesic orbit. Depending upon whether X is determined from the Jacobi 
or the generalized Jacobi equation, the results of this approach may be interpreted 



Tidal Dynamics in Kerr Spacetime 



4 



respectively in terms of the linear or nonlinear perturbations of geodesic orbits in the 
background gravitational field. 

Regarding the source of the gravitational field and the reference trajectory, we 
assume that m is on a stable equatorial circular geodesic orbit about a Kerr source. 
This choice for the orbit to some extent complements our previous work on the motion 
of astrophysical jets 0111, where the reference trajectory is a radial geodesic along the 
rotation axis of a Kerr source. The geodesies of Kerr spacetime have been discussed by a 
number of authors; see, for instance, ref. [Zj and references therein. The characteristics 
of the source and the reference trajectory will be described in detail in the next section. 
Various aspects of tidal dynamics in black-hole spacetimes have been discussed in a 
number of papers; see, for instance |H1 13 QHl ^1 ^] and the references cited therein. 
Furthermore, we treat free-particle orbits for the sake of simplicity; a corresponding 
discussion of black-hole accretion disks is beyond the scope of this investigation. 
Section El is devoted to the solution of the Jacobi equation. We show in section HI that 
a rotation introduced in section El leads to a complete transformation of the generalized 
Jacobi equation to an autonomous form. This equation is then analyzed in detail. 
Finally, section El contains a brief discussion of our results. 



2. Reference Worldline 



The spacetime region of interest is the exterior Kerr domain with the metric 

2 Mr 

+ —^{dt-asm^ed(f)y, (11) 

where a = J/M is the specific angular momentum of the source. Here [t, r, 6, 0) are the 
standard Boyer-Lindquist coordinates and 

E = + cos^ ^, A = r^ - 2Mr + . (12) 

The reference trajectory is assumed to be a stable circular orbit of fixed radius r in the 
equatorial plane 6 = tt/2; in terms of r, such orbits exist from infinity all the way down 
to the last stable circular orbits given by 

QM [M 

1 ±8aW^-3^ = 0. 13 

Throughout this paper, the upper (lower) sign refers to orbits where m rotates in the 
same (opposite) sense as the source. For r less than the solution of (O, there are 
unstable circular orbits that end at the null circular orbits given by 

3M [M 

— ±2a^7F=0- (14) 
It is useful to define the Keplerian frequency for the orbits under consideration here 

as 




ujk = ±J—. (15) 
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Let us note that in the orbital equations (fT^ and (HD), a prograde orbit becomes 
retrograde and vice versa when a —a; this circumstance explains why the combination 
auk usually appears in orbital equations. We assume that a > throughout; therefore, 
the sign of uk indicates the sense of the orbit. It follows from the geodesic equation 
that along the reference worldline 

i = ^^^, = ^^, (16) 

where 



N = ^1-^ + 2aujK (17) 

is such that the null circular orbits are solutions of = and we have assumed in ()16p 
that t = r = Oat0 = O. The constants of the motion for the reference worldline are 
the specific energy E and orbital angular momentum L, associated respectively with 
the timelike and azimuthal Killing vectors dt and d^, and are given by 

+ (18, 

To determine the orthonormal tetrad x, let us first consider the tetrad A.^, 
carried by the fundamental static observers in the exterior Kerr spacetime. In the 
equatorial plane and in terms of {t,r,6,(j)), the natural orthonormal tetrad of these 
observers along the coordinate directions is given by 



A% = 0, 0, j , A'^^,) = 1^0, 0, j , (19) 

f 1 \ .a ( 2Ma A\ 

A%=(0.0,;,0], A",, = (--^,0,0,-^), (20) 



where 



Next, we subject this tetrad to a Lorentz boost with speed /3, A^^.^^ — > A'^^^^, such that 
(0) ~ •^%) vxAX vector tangent to the worldline of the reference particle m. The 

new orthonormal tetrad along the orbit of m is 

A%)=7[A%+M%], Af,)=A%, (22) 
A% = A^(,), A'^(3)=7[A'^(3)+/3A%], (23) 

where 7 is the Lorentz factor corresponding to /5. We find from A'^j^q^ = A^j-q^ that the 
Lorentz pair (/3, 7) is given by 

/3=^ 7 = ^ (24) 
^ EN ' ^ A ^ ' 

Clearly, /? is positive (negative) for prograde (retrograde) orbits; moreover, 7 diverges at 
the null orbits (A^ = 0). The spatial triad of A'^j.^^ is along the spherical polar coordinate 
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Figure 1. Schematic drawing tliat depicts the paraUel transport of the spatial triad 
along a prograde circular orbit. 



directions, which are the radial, normal and tangential directions with respect to the 
orbit, and therefore needs to be rotated back such that the resulting tetrad would then 
be parallel propagated along the orbit as illustrated in Figure [H Thus 

^^"(2) = ^^(2) ' ^^(3) = ^^(1) sin V + A''(3) cos r]. (26) 
It follows from the vanishing of the covariant derivative of A^^^-^ along the orbit that 

V = ^kt, (27) 

where we have chosen the integration constant such that r/ = at r = 0. Note that 
(j) — rj = {N~^ — l)ri, which, to first order in M/r -C 1 and a/r <^ 1 can be written as 

<P-V^il^u^-^^)r- (28) 

that is, the difference between these angles is due to a combination of geodetic 
and gravitomagnetic precessions. The final result for the tetrad frame in [t, r, 6, (p) 
coordinates is then 

A^(o) = ^(l + aa;K,0,0,u;^), (29) 
A^j^^^ = --j= (^~^ sin r], — cos r], 0, —E sin rj^ , (30) 

fO'O'-'OV (31) 



r 
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(3) 



A 

L cos rj, — sin ?7, 0, cos rj 



(32) 



Finally, we need to compute the projection of the Riemann tensor along the orbit 
on this tetrad frame as in equation (@)). The result may be expressed as a 6 x 6 matrix 
(JZab), where A and B are indices that belong to the set {01,02,03,23,31,12}. This 
matrix has the form 



7^ 



S 

n 



n 

-£ 



(33) 



where £ and Ti are 3x3 symmetric and traceless matrices containing respectively the 
electric and magnetic components of the Riemann tensor. It turns out that in the case 
under consideration 



£ 



K 




k' 





k2 




k' 


ks 



H 



K 





h 





h 



h' 





h' 





where 



M 



K 



(34) 



(35) 



The proportionality 7?. oc /« is expected on the basis of the correspondence with 
Newtonian tides. In (jMj) . k2 = — (fci + ^3) is constant and is given by ^2 = — 2, 
while 



37^ cos^ rj, 



37^ sin^ rj, 



ki = l- 

h = COST], h' = —S'j^Psmr]. 
Here (/3, 7) is a new Lorentz pair, i.e. 7 



k' = —37^ sin rj cos rj, 



given by 



P 



r ujk — a 



A 



7 



rN 



(36) 
(37) 



(38) 



which reduces to (/3, 7) in the Schwarzschild limit (a = 0); moreover, 7 diverges at the 
null orbits (A^ = 0). 

Let us note that for the orbits under consideration and a < M, (3 G [—1/2, 1/2]; 
that is, (3 is positive (negative) for a prograde (retrograde) orbit and far from the source 
(r > 2M), \(3\ ~ ^M/r, so that /5 — >■ as r — > 00, then \f3\ monotonically increases 
from zero to 1/2 at the last stable circular orbits. Furthermore, for a > M, (3 is always 
negative for retrograde orbits, but it is not always positive for prograde orbits; in fact, 
j3 vanishes for a prograde orbit with radius r = /M. 

These curvature components can now be used in equations ((Zj) and to study the 
motion of the probe relative to m. It should be mentioned that our results for the tidal 
tensor £ are consistent with the work of Marck [Oj . 
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3. Jacobi Equation 

Let us now turn to the solution of tlie Jacobi equation 0. For simplicity, we use instead 
of the Fermi coordinates X" = (T, X), the dimensionless Fermi coordinates [f],^), 

cu^X":=(r7,x), (39) 

where x = (x, y, z). The Jacobi equation then takes the form 



drf 
drf 



X 


= S 


X 


s = 


z 




z 





37^ cos^ ?7 — 1 37^ sin rj cos r] 
37^ sin 7] cos 37^sin^r7 — 1 



0, 



cr 



(40) 
(41) 



Thus normal to the equatorial plane, we have a simple harmonic motion of frequency 
|a;i4:|(T, where a ranges from 1 to a/2 for a < M; that is, it increases monotonically from 
unity at r = 00 to at the last stable circular orbits. The equations of relative motion 
in the equatorial plane can be put into an autonomous form by a rotation to the radial 
and tangential directions (cf. Figure (H). That is, let 







X 




= R 


R = 


_c_ 




z 





cos rj sm rj 
— sin rj cos rj 



(42) 



Then, equation (jlUj) reduces to the autonomous system 

e-2C-37'e = 0, C + 2i = 0, (43) 

where an overdot denotes differentiation with respect to rj. Assuming that at 77 = 0, 
X = and the initial velocity of the probe (Vq ^ 1) has components 

Vo = Vo(sin cos (p, cos ^, sin sin (p) (44) 

in the local {X, Y, Z) system, we find 
Vo sin , 



y 
C 



A2 



■ [2(1 — cos A?7) sin + A sin Xrj cos (p] , 



a 



-Vo cos "i? sin arj, 

Vosin?? 2 4 , 

— — (37 rj — — sm Xrj) smy? + 2(1 — cos Xrj) cosy? , 

A^ A 



where A = a/4 — 37^. We note that A > 0, since 



A^ = 4^(1-6— + 8ac<;i^ 



(45) 
(46) 
(47) 

(48) 



m ■ r 

so that for a < M, X ranges from zero to unity; that is, it decreases monotonically from 
unity at r = 00 to zero at the last stable circular orbits. If the reference trajectory is 
one of the last stable circular orbits, then 7^ = 4/3, ?/ is given by (jlUj) with a = \pl and 

,^ = Vo sin '(9(?7 cos + r/^ sin (y9), (49) 
2 



C = K)Sint9[(r7 - -rf 



sin — ?7 cos Lp\ . 



(50) 
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Recalling the restrictions |X| ^ TZa and |V| -C 1, which in this case translate to |x| ^ 1 
and |x| <^ 1, we find that the validity of equations (07j), ()49|) and ()50|) is limited in time 
due to the presence of secular terms. 

It is interesting to note that terms that appear in the solution of the Jacobi equation 
in equations (jlHjl and (jTfjl with proper frequency X\ujk \ have frequency X\uK\N/{l+aujK) 
with respect to the coordinate time t; to lowest order in M/r <^ 1 and with a = 0, the 
deviation of this frequency from the Keplerian frequency corresponds to the Einstein 
precession frequency for an orbit of vanishing eccentricity, while the first-order correction 
in a/r ^ 1 corresponds to the de Sitter-Lense-Thirring precession frequency 

Finally, a remark is in order here regarding the boundary condition that X = at 
r = 0, which corresponds to our assumption that the probe is initially launched from 
m. In fact, it is only necessary that the probe be near the reference particle m, so that 
|X| must be initially very small compared to TZa] for the sake of simplicity, we choose 
|X(0)| = throughout this paper. This means that in terms of {^,y,C)y we always 
assume that 

(e,i/,c)u = o, (e,2/,c)u = Vo. (51) 



4. Generalized Jacobi Equation 





-2C 




= D 




c 


+2^ 


-c 




c 



The reduction of the Jacobi equation to the autonomous system (jl^ is basically due to 
the fact that RSR~^ = D, where D = diag (87^ — 1, —1) is a constant matrix. System 
()40|) then takes the form 



(52) 



which is equivalent to system ()43|) . The transformation from {x,z) to (^, C) involves 
a rotation from a local inertial frame to rotating axes; therefore, in equation we 
note the presence of Coriolis and centripetal terms on the left-hand side of this system. 
Moreover, this transformation is related to Hill's contributions to the classical three- 
body problem [121 HZ] ; indeed, system (jlH)) in the Newtonian limit 7 ^ 1 is equivalent 
to a limiting form of the Hill system discussed in ref . ^Hl , namely, system (PT|1 of ref . 
with k = ei = 0. 

It can be shown, after much algebra, that under the same rotation as in equation 
(jl21), the nonlinear generalized Jacobi equation can be rendered autonomous as 
well. This remarkable fact is due to the special symmetries of the Kerr spacetime 
Introducing = ^ — ( and = + ^, where 



A, 



R 



X 

z 



(53) 



we find that xx + zz = C,A.^ + CA(, xz 



zx 



eAc - and 



_ 2C - 372^ = 2A^[CA^ - {a' + 1)^A^] + 6^'^^ 
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-2y\^'f3CA^-^0, (54) 
y + a'{l - 2f)y = ^y[A| - {a' + l)Al - 9^'f3A^] - 2y[^'(3A^{^A^ - (A^) 

+ (37^ - ^)eA5 - (7^ + ^)CAd + 27W(2/Ac - Cy), (55) 
C + 2i = 2Ac[CAc - (a^ + - 6f /J^A^ - 2A^(f /^A^ + la'^A^ - (A^) 

+ l)y{yA^ - (y) + 2yy[f /3(A^ + 3) + a'A^] - 2^^Py\A^. (56) 

Starting with the exact solution of the Jacobi equation as the unperturbed solution, 
it is straightforward to develop a solution of equations (IM jl - lj^Uj) via the standard 
perturbation expansion in terms of Vq ^ 1. However, the results of such a perturbation 
scheme, based on the small parameter Vq, should be used only in conjunction with the 
higher-order tidal terms — that is, terms that have been neglected in the expansions 
of the gravitational potentials in equations (0)-® and would result in terms in the 
equation of motion (0) that go beyond the linear order in |X|. 

It follows from equations (HHI), jSHl) and ^ that the trajectory of the 

probe in Boyer-Lindquist coordinates, = itp,rp,9p,(l)p), is given by 

tp{T) = ^(1 + auK)T + C, (57) 



rp(r)=r + -l^e, (58) 
W = f + -^?/, (59) 
0p(r) = ^^i^r + ^^C. (60) 

Using 1^(77), y{rj) and C,{vi) from the (generalized) Jacobi equation and replacing r] = ujkT 
by ujkt, we can determine the path of the probe in the standard Kerr coordinate system. 
This is due to the fact that the generalized Jacobi equation is linear in the distance |X| 
away from the reference geodesic; therefore, X^A^j^^-j may be viewed as a generalized 
Jacobi field that is defined along the reference geodesic. It turns out that for small- 
amplitude perturbations (Vq ^ 1), the radial (^) and vertical {y) motions in general 
contain the basic proper epicyclic frequencies \\uik\ and cxltuxl, respectively. 

The autonomous system of equations naturally splits into equations for 

the vertical and equatorial motions. That is, if the probe is launched vertically in the y 
direction, then its motion will be confined to the y direction according to the generalized 
Jacobi equation. Similarly, if the probe is launched in the orbital (^, Q plane, then the 
probe will remain in this plane throughout its motion. 
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Vertical Motion 

If the probe is launched in the purely vertical direction relative to the reference orbit, 
then ^{rj) = ({rj) = for all 77 and hence the equations of motion reduce to 

y + a\l-2f)y = 0. (61) 

The vertical motion is simply uniform for y = ±l/\/2; that is, equation ()6H1 exhibits 
the critical speed Vc = 1/a/2 beyond which the character of the motion is opposite to 
the low-speed limit that agrees with Newtonian expectations. If the probe is launched 
with Vq < Vc, then it will begin to decelerate in agreement with Newtonian gravity. 
However, for Vq > Vc, the probe will accelerate, as treated in detail in [21 In fact, 
the vertical motion can be described exactly [21 HI E] , since equation (|HT|) implies that 

f = Vc' + {V,'-Vc')e''^'y\ (62) 

whose solution can be determined by quadrature. This solution represents timelike 
geodesic motion if the following condition is satisfied 

For Vq < Vc, condition (jUHj) is always satisfied and it follows from (jU^ that the motion 
of y is periodic and confined to the interval [—ymax, Z/max] , where 




On the other hand, for Vq > V^ the particle is accelerated to (almost) the local speed of 
light. In this way, the particle gains enormous tidal energy. According to condition (jUHj) . 
the range of y over which this happens monotonically shrinks to zero as Vq approaches 
unity. The question of the admissibility of Fermi coordinates in such a case is a difficult 
one; an unsuccessful attempt to answer such a question is contained in our recent 
work jH]. It is essential — for astrophysical applications — to connect these local results 
to what distant observers would measure [Sj. Indeed, the gravitational tidal energy of 
particles within the Fermi coordinate system can be transferred to the outside world 
(i.e. beyond the Fermi system) through collisions with other particles or emission of 
radiation. 

In our previous work involving tidal dynamics about a reference escape trajectory 
along the axis of rotation of the central source [Hj, we found that for relative speeds 
above the critical speed Vc, tidal deceleration occurs in a cone of angle G measured from 
the rotation axis, where G is given by tanG = l/K and corresponds to an angle of 
about 55°, while tidal acceleration occurs outside this cone, i.e. within the complement 
of the critical- velocity cone involving latitudes from — G' to G', where G' = 7r/2 — G 
and corresponds to an angle of about 35°. Based on these previous results, we expect 
that in the present situation the phenomenon of tidal acceleration for speeds above the 
critical speed may not be restricted to the vertical axis alone. To investigate this issue 
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in the case under consideration, we have integrated the full system of equations (jHD)- 
(j^mi with initial conditions (j^Tj) . Our results indicate that vertical acceleration — either 
parallel or antiparallel to the Kerr rotation axis — occurs when the probe is launched 
with Vq > I/V2 in a direction (-(9, ip), where for a given ip G [0, 27r), is restricted to a 
certain interval, namely, e [0, t^o)- To clarify the situation analytically, let us first note 
that if C) is a solution of equations then so is —y,()- In particular, 

equation (pT|) is clearly invariant under the parity transformation y \—>- —y; therefore, 
vertical motion is expected to occur symmetrically with respect to the equatorial plane. 
For this reason, we limit our considerations to y > 0, which corresponds to motion 
antiparallel to the Kerr axis in our convention (see Figure Q). Furthermore, let us note 
that with initial conditions (j^Tj) . vertical acceleration vanishes initially, i.e. y{r]) = at 
?7 = 0; therefore, it proves useful to compute y (0) using equation (j55|l . It turns out 
that 

y (0) = 2VoW{'d, ip) cos ^9, (65) 
where is a quadratic function of p = sin i} given by 

W = Co- VoCip - V^C2p'. (66) 
Here Cq, Ci and C2 are defined by 

C, = a\V^-]^), Ci = 37^/3 sin yp, (67) 
C2 = (2or2 + 1) _ (^2 + 2) sin^ ip. (68) 

For 1!} = 0, y (0) is positive for Vq > in fact, this is the case for an extended 

domain about the Kerr axis. The boundary of this domain is characterized hj W = 0, 
i.e. p{ip) = sint^o given by 

Pif) = ^(Vcn4QA - CO, (69) 

and p = 1 (i.e. = 7r/2), where y (0) vanishes, as illustrated in Figure El where 
Vo = 0.8 and (3 = 0.45. Let us recall that for a prograde reference orbit of radius r 
about a Kerr black hole, /3 increases monotonically from zero at r = 00 to 0.5 at the last 
stable circular orbit. For a fixed P G (0, 0.5] and Vq > Vc, the tidal-acceleration domain 
expands as Vq — > 1. On the other hand, for a fixed Vq G {Vc, 1), the acceleration domain 
in Figure 121 remains practically the same as (3 —>■ 1/2. For P —>■ 0, however, the shape 
of the domain changes significantly. In fact, for (3 = 0, it follows from equation ()69j) 
that p = po/l cosip\, where po defined by V3 po = y^l — V^ /Vq is a positive constant 
such that Po G (0, l/-\/6). Therefore, the acceleration domain in this case is within 
the intersection of the vertical strip defined by — po < x < po, where x = p cos ip> and 
y = psinyj, and the circle of unit radius centered at the origin of {x,y) coordinates. 

Though we have considered a simple circular equatorial orbit as the reference 
trajectory in this paper, it is not difficult to imagine that the tidal acceleration process 
described above — within the context of the physics of an accretion disk about the 
source — could possibly contribute to the formation of astrophysical jets. Moreover, 
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Figure 2. Plot of the boundary of the domain — in the {x, y) plane — in which the 
vertical tidal acceleration of the probe is positive for Vb > 1/ V^- The probe is launched 
in the direction ({}, (p) with respect to the local frame of the reference observer. We 
define x := p cob ip, y := pshup, where p = sin??. For this plot, Vb = 0.8 and [3 — 0.45. 
The upper part of the boundary curve is given by equation H69|l . The lower part of the 
boundary curve is an arc of the circle of radius p = 1 on the interval 'K+a < ip < 27r— a, 
where a ~ 0.94 radian corresponding to about 53°. 

tidal acceleration occurs symmetrically with respect to the equatorial plane of the Kerr 
source. This feature would be consistent with the observed occurrence of double jets, 
namely, a pair of relativistic outflows in opposite directions along the rotation axis of 
the central source. 

Equatorial Motion 

If y{ri) = for all rj, then the motion is confined to the (^, () plane and represents a 
timelike geodesic if 

^ = 1 - A| - - + i)e +e- ^I'pm^c - c\) - ^^'(^Ac - c^^f > o.(7o) 

A solution of the generalized Jacobi equation (GJE) in the (^, () plane is characterized 
by the parameters P, Vq and ip. For Vq extremely small compared to unity, the solution 
of the GJE is given essentially by the corresponding solution of the Jacobi equation as 
in the previous section. In general, as Vq slowly increases toward unity, the infiuence of 
the nonlinear terms in the GJE cannot be ignored. In fact, since ( in the Jacobi case 
contains a secular term, namely, —3{'~f/X)'^{Vosmi)sirnp)r] with i) = 7r/2, the solution 
of the GJE in general grows rapidly and leaves the admissible region of the Fermi 
coordinate system. This difficulty can be avoided, however, in the case of the purely 
radial variations with = or vr; in either case, the Jacobi equation has periodic 
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Figure 3. Plots of ^) versus 77, top panels, and Q versus t], bottom panels, for 
the generalized Jacobi equation with f3 — 1/5. The initial conditions are such that the 
motion is confined to the equatorial plane with a radially inward initial velocity such 
that i9 — tt/2, If ^ tt and Vq = 0.01. Note that the fast oscillation in the above panels 
has essentially the same frequency as in the corresponding solutions of the Jacobi 
equation, (^j, O), given by O = -(Vo/A) sinA?7 and O = (2Vo/A^)(l - cos At?). 



solutions corresponding to the stability of the circular reference orbit under small radial 
perturbations. The solution of the GJE may then contain quasi-periodic oscillations 
exhibiting a complex beat phenomenon involving several frequencies as demonstrated 
in Figure El However, such oscillations appear to occur only for Vq -C 1, in which 
case equations are generally not adequate physically, since higher-order tidal 

terms of comparable magnitude have been neglected in our analysis. Even if a quasi- 
periodic oscillation of the type illustrated in Figure |H1 survives the inclusion of relevant 
higher-order tidal terms, the amplitude of the effect may not be large enough to be 
of any physical significance in connection with the observed quasi-periodic oscillations 

(QPOs) y^. 

5. Discussion 

We have studied the generalized Jacobi equation for a circular reference geodesic orbit in 
the equatorial plane of the exterior Kerr spacetime. This equation has been reduced to 
an autonomous system; the corresponding tidal dynamics can then be naturally divided 
into vertical and equatorial motions relative to the circular orbit. In connection with 
vertical motion, we have clarified the role of the critical speed Vc = 1/ in the character 
of the solutions of the generalized Jacobi equation. 

The general approach developed in this work may be of interest in connection with 
relativistic physics inside a space station, the tracking of artificial Earth satellites as well 
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as satellite-to-satellite Doppler tracking. Moreover, in connection with the relativistic 
astrophysics of accretion disks around (rotating) astronomical sources, our results may 
be relevant for the gravitational aspects of the complex plasma physics that would be 
involved in the formation of high-energy jets [20] • That is, the vertical acceleration 
phenomenon explored in this work may be related to how relativistic jets get started 
above and below an accretion disk around a Kerr black hole. Once the flows are 
magnetically confined to regions near the axis of rotation of the central source, the 
results of previous investigations El 112] involving tidal dynamics of ultrarelativistic 
particles would have to be taken into account for the parallel and antiparallel flows along 
the Kerr rotation axis. 

References 

[1] Synge J L 1960 Relativity: The General Theory (Amsterdam: North-Holland) 

[2] Chicone C and Mashhoon B 2002 Class. Quantum Grav. 19 4231 

[3] Chicone C and Mashhoon B 2005 Preprint 'gr-qc /0511129 

[4] Mashhoon B 2005 Int. J. Mod. Phys. D 14 2025 

[5] Chicone C and Mashhoon B 2005 Ann. Phys. (Leipzig) 14 290 

Chicone C and Mashhoon B 2004 Class. Quantum Grav. 21 L139 

Chicone C and Mashhoon B 2005 Class. Quantum Grav. 22 195 

Chicone C and Mashhoon B 2005 Astron. & Astrophys. 437 L39 

Chicone C and Mashhoon B 2005 Ann. Phys. (Leipzig) 14 751 
[6] Chicone C, Mashhoon B and Pimsly B (2004) Int. J. Mod. Phys. D 13 945 

Chicone C, Mashhoon B and Punsly B 2005 Phys. Lett. A 343 1 
[7] Chandrasekhar S 1992 The Mathematical Theory of Black Holes (New York: Oxford University 

Press) 

[8] Mashhoon B and Theiss D S 1982 Phys. Rev. Lett. 49 1542 

[9] Marck J-A 1983 Proc. Roy. Soc. London A 385 431 
[10] Bini D, Gerahco A and Jantzen R T 2005 Class. Quantum Grav. 22 4729 
[11] Ishii M, Shibata M and Mino Y 2005 Phys. Rev. D 71 044017 
[12] Kojima Y and Takami K 2006 Glass. Quantum Grav. 23 609 

[13] Kato S, Fukue J and Mineshige S 1998 Black-Hole Accretion Disks (Kyoto: Kyoto University 
Press) 

[14] Ahev A N, Gal'tsov D V and Petukhov V I 1986 Astrophys. & Space Set. 124 137 
[15] Ahev A N and Gal'tsov D V 1988 Sov. Astron. Lett. 14 48 
[16] Chicone C, Mashhoon B and Retzloff D G 1996 Ann. Inst. H. Poincare 64 87 
[17] Chicone C, Mashhoon B and Retzloff D G 1999 Helv. Phys. Acta 72 123 

[18] Carter B 1973 Black Holes edited by C DeWitt and B S DeWitt (New York: Gordon and Breach) 
[19] Strohmayer T E 2001 Astrophys. J. 552 L49 

Strohmayer T E 2001 Astrophys. J. 554 L169 

Wagoner R V, Silbergleit A S and Ortega-Rodriguez M 2001 Astrophys. J. 559 L25 
[20] Punsly B 2001 Black Hole Gravitohydromagnetics (New York: Springer- Verlag) 



